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RIEMANN-HILBERT PROBLEM AND 
THE DISCRETE BESSEL KERNEL 

Alexei Borodin 



Abstract. We use discrete analogs of Riemann-Hilbert problem's methods to de- 
| rive the discrete Bessel kernel which describes the poissonized Plancherel measures 

for symmetric groups. To do this we define discrete analogs of a Riemann-Hilbert 
problem and of an integrable integral operator and show that computing the resolvent 
of a discrete integrable operator can be reduced to solving a corresponding discrete 
Riemann-Hilbert problem. 

We also give an example, explicitly solvable in terms of classical special functions, 



when a discrete Riemann-Hilbert problem converges in a certain scaling limit to a 
conventional one; the example originates from the representation theory of the infinite 
symmetric group. 



(N 
> 

', 1. Introduction 

O ■ 

Let S(n) be the symmetric group of degree n. It is well known that its irreducible 
representations can be naturally parametrized by all the partitions A = (Ai > A2 > 
• • • > Afc) of the number n into natural summands: n = Ai + A2 + • • • + A&, see, 
e.g., [JK]. We shall denote the set of all partitions of n by V n . For every partition 
A, let dim A be the dimension of the irreducible representation corresponding to A. 
The well-known Burnside formula for finite groups implies that 



X 



dim 2 A 



n\. 



Thus, we can construct a probability distribution M n on V n by setting the weight 
of A G V n equal to 



M n (X) 



dim 2 A 



nl 

The M n 's are called the Plancherel distributions. These distributions can be ob- 
tained from the uniform distributions on symmetric groups via the Robinson- 
Schensted correspondence, see, e.g., [Sch]. In particular, the distribution of the 
largest part Ai of the random partition of n with respect to the Plancherel distribu- 
tion coincides with that of the longest increasing subsequence of random uniformly 
distributed permutation of degree n, see, e.g., [BDJ1, Appendix]. 

Let us organize a new probability distribution M 6 on the set V = U n> iV n of all 
partitions depending on a positive parameter 9 > as follows. Let A G V n . Set 

on 

M e (\) = —e- e -M n (\). 
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As was shown in [BOO], the correlation functions of a point process naturally 
attached to M e have determinantal form with a certain kernel expressed through the 
J Bessel functions. This powerful fact allowed us to study asymptotic properties of 
the Plancherel distributions and, in particular, to prove the Baik-Deift-Johansson 
conjecture [BDJ1], [BDJ2], about the asymptotic behavior of finitely many first 
parts of the large random partition. 1 The same correlation kernel has also arisen in 
[J] in the asymptotics of orthogonal polynomial ensembles related to the Plancherel 
measures. See also [B013] for a discussion of connections between the approaches 
of [BOO] and [J]. 

The appearance of Bessel functions in [BOO] seemed rather mysterious. A nice 
representation theoretic explanation was suggested in [Ok] — for a much wider class 
of measures on partitions it was shown that the correlation functions are given by 
determinantal formulas with a kernel which has a certain integral representation. 
In the particular case of Plancherel distributions this approach leads to the integral 
representation of Bessel functions, see [BOk, §4]. 

In this note we will obtain the Bessel functions in another way. We will show that 
the correlation kernel can be defined as a solution of a certain "discrete Riemann- 
Hilbert problem" . 

First, we will define discrete analogs of Riemann-Hilbert problems and integrable 
integral operators and show how the resolvent of a discrete integrable operator can 
be obtained from a solution of the corresponding discrete Riemann-Hilbert problem. 
The result is parallel to the known result in the continuous case, see [IIKS] , [D2] . 

As was pointed out by P. Deift and A. Its, our setting of the discrete Riemann- 
Hilbert problems is similar to the pure soliton case in the inverse scattering method, 
see [BC], [BDT], [NMPZ, Ch. III]. 

Then, applying discrete analogs of standard methods of the conventional Rie- 
mann-Hilbert problem's techniques, we will obtain a system of differential equations 
on the matrix elements of the solution of our concrete discrete Riemann-Hilbert 
problem which will lead to the Bessel equation. 

We will also demonstrate an explicitly solvable example of a discrete Riemann- 
Hilbert problem (more general then the one that arises from Plancherel measures) 
and a continuous Riemann-Hilbert problem such that the discrete problem con- 
verges to the continuous one in a certain scaling limit. The solutions of these 
problems are expressed through classical special functions. The example has orig- 
inated from a problem of harmonic analysis on the infinite symmetric group, see 
[KOV], [BOll], [B012], [B013]. 

The paper is organized as follows. §1 is the introduction. In §2 we describe a 
problem whose solution provides the correlation kernel for the poissonized Plancherel 
measures (the discrete Bessel kernel) . In §3 we review general facts about integrable 
integral operators and Riemann-Hilbert problems. In §4 we define discrete analogs 
of Riemann-Hilbert problems and integrable integral operators and show how com- 
puting the resolvent of a discrete integrable operator is reduced to solving a discrete 
Riemann-Hilbert problem. In §5 we apply the general approach discussed in §3 to a 
special case which is relevant for us. In §6 we do the same in the discrete situation. 
In §7 we solve the discrete Riemann-Hilbert problem attached to the problem of 
§2 using discrete analogs of standard methods of continuous problems. In §8 we 



1 We refer to [BDJ1] and [AD] for references to earlier works on the subject. Note, in particular 

j 1 _ _ _______ .T n . r-v -I . rr/-nl 
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discuss the example of explicitly solvable discrete and continuous Riemann-Hilbert 
problems mentioned above. 

In the preprint version of this paper there was no general setting for discrete 
Riemann-Hilbert problems and discrete integrable operators. Only the special case 
of §6 was worked out. After the preprint had appeared, Percy Deift suggested a 
general approach to the discrete situation which we present in §4. I am very grateful 
to him for allowing me to reproduce his results in this text and for a number of 
valuable suggestions. 

I would like to thank Grigori Olshanski. Without his constant support and 
stimulating discussions, this work would never be done. I would also like to thank 
Alexander Its for explaining the basics of Riemann-Hilbert problems to me and for 
helpful comments. 

2. Setting of the problem 

We refer the reader to [BOO] for a detailed exposition of the material of this 
section. 

Let us associate to any Young diagram A a finite point configuration X(X) in U 
as follows. Denote by (pi, • • • ,Pd I 9i> • • • > 9d) the Frobenius coordinates of A and set 

v/n [1 11 1 

X(X) = <-qi - -q d - -,p d + J* •••>?! + g 

We define the correlation functions pk(%i, • • • , Xk), k = 1, 2, . . . , of the measure 
M e by 

Pfc(xi,...,x fc ) = M e {X | {x 1 ,...,x k } C X(X)}, x 1 ,...,x k e Z'. 

As was shown in [BOO], the correlation functions have determinantal form, 

Pk{x\, . . . , Xk) = det[K(xi, Xj)]i tj=1 , 

with a certain kernel K(x, y) on Z'. 

This kernel can be defined by the formula K = L(l + , where the kernel L 
(also on U) has a particularly simple form, see [BOO, Proposition 2.3]. Namely, its 
block form corresponding to the partition of Z' into the set of positive half-integers 
Z' + and negative half-integers Z'_ is 



L(x,y) = 



n e~ L_ 

r(x+i)r(-y+i) x-y 

— x + y 

e~^~ i n 

Lr(-x+i)r(y+i) x-y 



(2.1) 



The notation means that L(x, y) = if x and y are of the same sign; 



9^ 1 

L(x, y) = — — —y- , x G Z' + , y G Z'_; 

T(x+ ±)T(-y+ ±) x-y 



-x-\-y 

e^ 1 i 



L(x,y) = — r- — ; r- , x G Z' , y G Z', . 

1 ' yJ Y{-x + \)Y{y + \)x-y' + 

This fact is quite elementary. It is more difficult to obtain an explicit formula 
for K. In [BOO] we gave several such formulas. Here is one of them (Theorem 1 
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Theorem 2.1. Define the kernel K(x,y) on U = U + U U_ by the formula 



K(x,y) = Ve 



J i J . i —J . i J i J i J i + J . i J ,i 



x-y x-y 
l J i +J i J , i J i J l — </ l J i l 

x-y x-y 



(2.2) 



where J x = J x {2\/0) is the Bessel function of order x and argument 2\^9, and the 
diagonal entries are determined by the L 'Hospital rule. Then 

T, 

K = 



1 + 

where L is defined by (2.1). 

We call K(x,y) the discrete Bessel kernel. 

In [BOO] we gave no explanation why the Bessel functions appear in the picture, 
there we have just verified that the relation K = L(l + L) -1 holds. Several such 
explanations exist by now. Originally, we obtained the discrete Bessel kernel as a 
limit of the hyper geometric kernel introduced in [B012], see also §8. The hyperge- 
ometric kernel describes a certain two-parametric family of measures on partitions 
called z-measures which can be degenerated to the Plancherel measures. K. Johans- 
son independently obtained the restriction of the discrete Bessel kernel to Z' + as a 
limit of Christoffel-Darboux kernels for the Charlier polynomials [J] . A. Okounkov 
showed that the Plancherel measures as well as the z-measures fit into a much 
more general infinite-parametric family of measures on partitions, and he gave an 
integral representation for the corresponding correlation kernels [Ok]. In the par- 
ticular case of the Plancherel measures one obtains the integral representation of 
the Bessel function in this way, see [BOk, §4, Ex. 1]. 

In this note we suggest yet another way to obtain the discrete Bessel kernel. We 
will consider L and K as discrete analogs of integrable operators in the sense of 
[IIKS] , see also [D2] . To my best knowledge, such an object has not been discussed 
in the literature before. 2 

It is known that if L is an integrable operator in the usual sense then the operator 
L(l + L) _1 is also integrable, and it can be expressed through a solution of a 
certain Riemann-Hilbert problem (RHP, for short) associated with L, see [IIKS], 
[D2], and §3. We will define a discrete analog of a Riemann-Hilbert problem and 
we will show that solving a certain discrete Riemann-Hilbert problem (DRHP, for 
short) is equivalent to computing the resolvent of a discrete integrable operator. 
Furthermore, we will show how to obtain the discrete Bessel kernel by employing 
discrete analogs of standard methods used in handling conventional (continuous) 
Riemann-Hilbert problems. 



3. Integrable operators and Riemann-Hilbert 
problems. General approach 

This section gives a brief review of the formalism of integrable operators and 
corresponding Riemann-Hilbert problems. We shall follow [D2] in our description 
of the material. 

2 Of course, the Christoffel-Darboux kernels for discrete orthogonal polynomials can be consid- 
ered as discrete integrable operators, but, as far as I know, this observation has never been used 

_ _ r 
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Let E be an oriented contour in C. We call an operator L acting in L 2 (E, \d(\) 
integrable if its kernel has the form 

for some functions /j, gj, j = 1, . . . , N. We shall always assume that 

JV 

£/i(C)fc(C) = o, (gs, 

3 = 1 

so that the kernel L(£, £') is nonsingular (this assumption is not necessary for the 
general theory). 

We did not impose any restrictions on the functions fi,gi and on the contour 
E here, these restrictions depend on a particular problem. For example, one can 
demand that j^g% £ £ 2 (E, \d(\) fl L°°(E, \d(\), and the contour E is such that the 
Cauchy Principal Value operator H, 

{C'eS:|C'-CI>e> 

is L 2 -bounded. These restrictions imply, in particular, that the operator L is a 
bounded operator in L 2 (E, \dC\), see [D2]. 

The concept of an integrable operator was first introduced in [IIKS]. 

It turns out that for an integrable operator L such that (1 + exists the 

operator K = L(l + L)~ x is also integrable. 

Proposition 3.1 [IIKS]. Let L be an integrable operator as described above and 
K = L(l + L) _1 . Then the kernel K((, (') has the form 

K(CX) = E ' =1 ^p {C '\ CC'GE, 

where 

Fj = (1 + L)~ 1 fj, G j = (l + L t )- 1 g j , j = l,...,N. 

If Ef=i fAOQiiO = on E, then Ef =1 ^(0^(0 = on E as well. 

A remarkable fact is that Fj and Gj can be expressed through a solution of an 
associated Riemann-Hilbert problem. 

As we move along E in the positive direction, we agree that the (+)-side (re- 
spectively, (— )-side) lies to the left (respectively, right). 

Let v be a map from E to Mat(/c, C), k is a fixed integer. 

We shall say that a matrix function m : C \ E — > Mat(/c, C) is a solution of the 
RHP (E, v) if the following conditions are satisfied 



m(() is analytic in C \ E, 



• m+(C) = m_(C)v(C), C e S, where m±(C) = hm m(C), 

c'-c 

C'€(±)-side 



m(C) — >■ / as C — > oo. 
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Proposition 3.2 [IIKS]. Let L be an integrable operator as described above and 
m(() be a solution of the RHP (E, v) with 

v(0=I-27rif(C)g(C) t eM a t(N,C), 



where 



/ = (h, ■ ■ ■ , InY , g = (#!,..., #jv)*- 



Then the kernel of the operator K = L(l + L) 1 has the form 

K(U) = ^f^, C,C'€E, 

where 

F = (Fi, . . . , Fjv) , G = (G\, . . . , Gn) 

are given by 

F(() = m±(() /(C), G(C) = (m^C))" 1 g(Q. 

4. Discrete integrable operators and discrete 
rlemann— hllbert problems. general approach 

In this section we prove discrete analogs of the two results stated in the previous 
section. The material of this section is due to Percy Deift. 

Let X be a discrete locally finite subset of C. We call an operator L acting in 
£ 2 (X) integrable if its matrix has the form 

f T,f=ifj(x)gj(x') , 
L(x,x') = l ' x ^ x ' (4.1) 

for some functions /-,-, gj on X, j = 1, . . . , N, satisfying the relation 

N 

'£f j (x)g j (x) = 0, *e£. (4-2) 
i=i 

We will assume that fj, gj G £ 2 (X) for all j. 

We will also call an operator in £ 2 (X) integrable if its matrix has the form (4.1) 
for x 7^ x' and has arbitrary (not necessarily zero) diagonal values. 

Set 

/ = (/i,-..,/jv)*, 9 = (0i,...,&\r)*. 
Then (4.2) can be rewritten as g t (x)f(x) = 0. We will also assume that the operator 

( Tft >M = E ^ < 43 > 

is a bounded operator in £ 2 (X) 3 . These restrictions guarantee that L is a bounded 
operator in £ 2 (X); a proof of this fact is contained in the proof of Proposition 4.3 
below. 

As in the continuous case, it turns out that if the operator (1 + L) is invertible 
then K = L(l + L)~ x = 1 — (l + L) -1 is also an integrable operator with possibly 
nonzero diagonal entries. 
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Proposition 4.1. Let L be an integrable operator as described above and K = 
L(l + L)~ x . Then the matrix K(x, x') has the form 

K(x,x) = — , x^x, 

x — X 

where 

Fj = (1 + L)~ 1 f j , G j = (l + L t )- 1 g j , j = l,...,N, 

and 

K(x,x) = - L {x,y)K(y,x). 

yEX, y^x 

Moreover, X^/Li Fj( x )Gj(x) = for any x G X. 

Proof We will follow the lines of the proof of Proposition 3.1 from [KBI, XIV. 1] 
The defining relation K + LK = L reads 

^—^ x — t x — y 

Assume x ^ y. Multiplying both sides by (x — y) and using the relation 

x-y =1 t - y 
x — t x — t 

we get 

(x - y)K(x, y) + J2 f ^TT {t ~ y)K{t ' V) = ft{x) \ 9(y)-J29(t) K ^y) ] • 

t=£x V t=/=x 

Thanks to (4.2), the restriction t ^ x in the second summation above can be 
removed, and we obtain 

E -*) + L ^ *)) (* - y) K ^ y) = /'(*) (C 1 - Rt )y) (v)- 
t 

Applying (1 + L)~ x to both sides we see that (x — y)K(x, y) = F t (x)G(y), where 

F = (Fi, . . . , Fat)*, G = (G\, . . . , GnY- 

This proves the first claim of the proposition. 

Since the diagonal entries of L are zeros, the relation K + LK = L implies 

K(x, x) = - E L(x, y)K(y, x) + L(x, x) = - ^ L(x, y)K(y, x), 

y^x y^x 

and this proves the second claim of the proposition. 

The relation F t (x)G(x) = will be proved later, see (4.12). □ 
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Remark 4.2. It is not difficult to see that for an integrable operator V with 
arbitrary diagonal entries bounded from —1, and such that (1 + V) is invertible, 
the operator W = Vil + V) -1 is integrable. Indeed, if we denote by Vd the diagonal 
part of the matrix V and by Vq the off-diagonal part of V (that is, V = Vd + Vq) 
then it is easily seen that 

W = Vd(l + Vd)- 1 +W(l + V d )-\ 

where 

w= (i + v d )-iv 



(l + Ci + Vd)- 1 ^)- 1 

is integrable by Proposition 4.1. This implies that W is also integrable. 

Similarly to the continuous case, Fj and Gj from Proposition 4.1 can be expressed 
through a solution of an associated discrete Riemann-Hilbert problem. 
Let w be a map from X to Mat(/c, C), k is a fixed integer. 

We shall say that a matrix function m : C \ X — > Mat(/c, C) with simple poles at 
the points x G X is a solution of the DRHP (X, it?) if the following conditions are 
satisfied 

• m(C) is analytic in C \ X, 

• Res£ =x m(C) = lim (m(()w(x)) , iGl, 

• m(C) — > / as C — > oo. 

We will also call the jnm^> matrix. 

If the set X is infinite, the last condition must be made more precise. Indeed, 
a function with poles accumulating at infinity cannot have asymptotics at infinity. 
One way to precise the condition is to require the uniform asymptotics on a sequence 
of expanding contours, for example, on a sequence of circles \(\ = a^, — > +oo. 

In order to guarantee the uniqueness of solutions of the DRHPs considered below 
we will always assume that there exists a sequence of expanding contours such that 
the distance from these contours to the set X is bounded from zero, and we will 
require a solution m(Q to have the proper asymptotic behavior on these contours. 

The setting of the DRHP above is very similar to the pure soliton case in the 
inverse scattering method, see [BC], [BDT], [NMPZ, Ch. III]. 

Proposition 4.3. Let L be an integrable operator as described above and m(() be 
a solution of the DRHP (X, w) with 

w (x) = -f(x)g(xY G Mat (A/ - , C). 

Then the matrix K = L(l + L)~ x has the form 

G\x')F(x) 

K(x,x') = 



li 



x 7^ x', 



G t (x) lim (m'(() f(x)) , x = x', 



where m'(() = — rf^- , and 
d( 



F(x) = lim (m(C) /(*)) , G{x) = lim ((m^C))" 1 g(x)) 
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Lemma 4.4. If m(C) is a solution of the DRHP (X, w), and w 2 (x) = for some 
x G X then the function 

™(0(/ + f^) (4.4) 
is analytic in a neighborhood of x. 

Proof. In a neighborhood of x we have (777(C) has a simple pole at x) 

m(() = -^- + B x + 0((-x) 7 
Q-x 

where A x and B x are constant matrices. Then 

Res c=x m(C) = A x , m{()w{x) = + B x w(x) + 0(( - x). 

The residue condition of the DRHP implies that 

A x w(x) = 0, A x = — lim m(Qw(x) = —B x w(x). 



Hence, 



m (() = B x [-^+I)+0((-x) 



and, recalling that w (x) = 0, we get 

m(C) (l+^l) =B X + 0(1) 

which is analytic near x. □ 

Note that the jump matrix w(x) = —f{x)g{x) t satisfies the condition w 2 (x) = 
at every point x G X because of (4.2). 

Lemma 4.5. Let mi(() and 777.2(C) be solutions of DRHPs with the same jump 
matrix w(x), w 2 (x) = for all x G X, and possibly different asymptotic conditions 
at infinity. If 7772(C) ^ s invertible for C G C \ X then 7771(C) (7772(C)) -1 ^ s analytic in 
C. 

Proof. Since w 2 (x) = 0, the determinant of the matrix + ^z^r) ^ s identically 
equal to 1. In particular, this matrix is invertible. 

If 777i(C) and 7772(C) are solutions of DRHPs with the same jump matrix w(x) 
and 7772(C) is invertible then 

m 1 (C)(m 2 (C))- 1 =m 1 (C)(/+^) ' MO)" 1 

mi(C) [I+T^] m 2 (C) U + 



C-xJJ V. V C-x 

is analytic near x G X by Lemma 4.4. Since s 6 X is arbitrary, mi(C) (7772(C)) -1 is 
analytic in C. □ 
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Corollary 4.6. A DRHP with the jump matrix of the form w(x), w 2 (x) = for 
all x G X, and arbitrary invertible asymptotics at infinity whose any solution is 
invertible in C \ X has at most one solution. 

Proof. Indeed, if mi(C) and m2(C) are two solutions then ?77i(C) (7772(C)) -1 has no 
singularities by Lemma 4.5 and is asymptotically equal to I at infinity By the 
Liouville theorem, mi(C) = 7772(C)- D 

Lemma 4.7. Ifm(() is a solution of a DRHP with the jump matrix w(x), w 2 (x) = 
for all j£l, and asymptotics I at infinity then det m(C) = 1. Hence, such DRHP 
has at most one solution. 

Proof. Lemma 4.4 implies that 

w(x) 



det m(C) = det 



rn(C) I + 



C-x 



is analytic near x G X. Since s 6 lis arbitrary, det 777(C) is analytic in C, and the 
asymptotic condition 777(C) ~~ * I as C - * 00, implies that det 777(C) = 1- ^ 

Proof of Proposition 4-3. The proof is similar to that of Proposition 3.1 from [D2]. 
It is based on the following commutation formula: if P>i and B 2 are Banach spaces 
then for any bounded operators D : E>i —> B 2 and E : B 2 — > -Bi, 

A • I) tJ-^ 1 



DE+X ED+X 

in the sense that if (—A) 7^ lies in the resolvent set of ED then (—A) lies in the 
resolvent set of DE and (DE + A)" 1 = A -1 (l - D(ED + A)" 1 ^), see e.g. [Sak]. 
This simple result turns out to have a large number of applications in mathematical 
physics, see [Dl]. 

Let Rf denote the map of right multiplication by the column iV-vector / taking 
row iV-vector functions to scalar functions: 

N 

(Rfh)(x) = h{x)f{x) = ^ h i{x)fi(x), h = (h u . . . , h N ), 

i=i 

and let R g t denote the map of right multiplication by the row iV-vector g l taking 
scalar functions to row iV-vector functions: 

(R g th)(x) = h(x)g t (x) = (h(x)g 1 (x), . . . , h(x)g N (x)). 

By applying Rf to a iV x N matrix valued function we will mean the application of 
Rf to every row of the matrix and getting a column iV-vector function as the result. 
Similarly, by applying R g t to a column iV-vector function means the application of 
Rgt to every coordinate of the vector and getting aniVxiV matrix valued function 
as the result. 

Now observe that the operator L is of the form L = DE, where D = Rf, 
E = TRgt , and T was introduced in (4.3) . The operator D maps the space B\ of row 
iV-vector functions on X with coordinates from £ 2 (X) to B 2 = £ 2 (X), and E maps 
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This implies that L is a bounded operator in £ 2 (X) if the functions fj,gj are 
bounded and the operator T is bounded. Note that for the boundedness of L 
instead of the boundedness of T we could require the boundedness of E = TR g t . 

Let us apply the commutation formula for A = 1. Since (1 + L) is invertible, we 
have (using Proposition 4.1) 

F = (1 + L)- 1 / = (1 + DE)~ x f = f-D Ef 

= f-Rf TR a tf = f -R f Tfg* 

J ED + 1 9 J J ED + 1 J y 

= f -Rf EDI = f -R f ( 1 J I 

J ED + 1 J J V ED + lJ 

= RfiED + iy 1 !. 

Here J denotes the n x n matrix valued function which is identically equal to 1. 

Set fx = (ED + Note that all matrix elements of (// — I) are in £ 2 (X). We 

have just showed that 

F = (1 + L)- 1 / = Rf(ED + l)" 1 / = ^, (4.5) 

where \i = I — ED fj,, or 

/*(*) = !- £ ^KMM. (4.6) 

x y 

yEX,y^x 

Since 

can be obtained from L(x,x') defined by (4.1) by the change / — > —g, g — > /, we 
get 

G=(l + L t )- 1 ^ = /I^ (4.7) 
where _ 

^)=/+ ^ m(s f t){y) . (4.8) 



Set 



x — y 

VEX, y^x 



m(Q=I-T, ^W M ^ m(() = J + ^ MgM . (4.9) 

cf. (4.6), (4.8). We have 

Res c=x m(C) = (/(/*)(*) = - Jim MO (/</*)(*)) . 

Indeed, 
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as (fg^ix) = 0by (4.2). 

Thus, m(C) is a (unique by Lemma 4.7) solution of the DRHP (£, w) with 
w(x) = —f(x)g(xY (the asymptotic condition immediately follows from (4.9)). 
Similarly, using (4.5) we get 

F(x) = fi(x)f(x) = lim (m(C)f(x)) . (4.10) 

Same arguments show that fh(() is a solution of the DRHP (X, g(x)f(x) t ) : and 

G(x) = Ji(x)g(x) = lim (m(£)gr(a;)) . (4-11) 

By Lemma 4.4, the functions 

m(0(/ + /( y) and ™(f) (/- 

are analytic near x. Then, using the relation (f(x)g(x) t ) 2 = 0, we have 

m(0««(0 = m(0 (/ + ll f^f) (i - ^f) «'(0 

- (»<o(^W)) W))' 

which is analytic near x. Since x 6 I is arbitrary, m(C)m t (C) is analytic in C, and 
the asymptotic condition of DRHP implies that m*(C) = (m(£)) _1 . Thus, (4.11) 
turns into 

G(x) = lim ((m*(C))-^(x)) . 

Q^x 

Note that by (4.10) and (4.11) we have 

G(xYF(x) = lim {gixYimiOr'miOfix)) = g(xYf(x) = (4.12) 

which proves the last claim of Proposition 4.1. 

Finally, on the diagonal by Proposition 4.1, (4.2), (4.5), (4.9) we get 

*(*,*) = - E H*,v)K<y,*) = - E 9{v) ' ,ix) G(X) ' F{V) 

^-^ x — y y — x 

y€X,y^x yEX,y^x 

=<*»)'( E f^)/w=GW'( E £ ^ i f^V> 

y/x / \j/€£, y/x / 

= G(xY lim(m'(C)/(*)). □ 
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5. INTEGRABLE OPERATORS AND 
RlEMANN— HlLBERT PROBLEMS. SPECIAL CASE 

In this section we apply the general formalism of §3 to a much more special 
situation. 

Let E = £/ U S/7 be a union of two contours, and assume that the operator L 
in the block form corresponding to this splitting is as follows 



L(x,y) 



q h I (x)h II {y) 
x-y 

h I (y)h II (x) q 

x-y 



(5.1) 



for some functions hj( ■ ) and hn( ■ ) defined on E/ and £//, respectively. 
Then the operator L is integrable with N = 2. Indeed, 

L(x,y) = fl{x)9l(y) + Mx)92(y) , z,yGE, 
x-y 

where 

r x e Ej, r 0, x e Ej, 

[0, £ G E/j, [ hn(x), x G E/j, 

( 0, x G £/, f /i/(x), x G E 7 , 

[ hn{x), x G E/7, [0, x G E//. 

Then the jump matrix of the corresponding RHP has the form 

1 -2nih 2 j(x) 



x G £/, 

'•(•'•) = < , , ' (5-2) 

-27Tih 2 n (x) J' XGE/7 ' 

It can be easily seen that the RHP in such a situation is equivalent to the following 
set of conditions: 

• matrix elements mn and TO21 are holomorphic in C \ E77; 

• matrix elements ui\i and m22 are holomorphic in C \ Ej; 

• on E/j the following relations hold 

mn + (x) — mn_ [x] = —2m hj(x)mi2(x), 
m 2 i + (x) - m 2 i-(x) = -2ni h 2 I (x)m 2 2{x)] 

• on E/ the following relations hold 

mi2 + (x) — mi2-(x) = —2-ni hj I (x)mn(x), 
m 2 2 + {x) - m22-{x) = -2ni hj I (x)m2i(x); 

• m(x) ~ 1 as x — > 00. 

According to Proposition 3.2, the kernel K(x, y) in the block form corresponding 
to the splitting E = E/ U E// looks as follows 

K{x,y) = 

h I {x)h I {y){-m 11 {x)m2i{y)+m2i{x)m 11 {y)) h j {x)h I j (y) (mn {x)m 2 2 (y) — m 2 i (x)m 12 (y)) 

x—y x—y 

h II (x)h I (y)(m22(x)mii(y)-m 12 (x)m2i(y)) hg (x)hu (y) (-m 22 (x)m 12 (y)+m 12 (x)m 2 2 (y)) 
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6. Discrete integrable operators and discrete 
rlemann— hllbert problems. special case 

Similarly to §5, we apply the general approach of §4 to a special case. 

Let X be a discrete locally finite subset of C and let X = Xi U Xu be its splitting 
into two disjoint parts. Let hi(-), hu{ ■ ) be two functions defined on Xj and Xjj, 
respectively. We will assume that hi G l 2 (Xi), hu G £ 2 (Xu). 

Consider a matrix L of size X x X which in the block form corresponding to the 
splitting X = Xi U Xu looks as follows, cf. (5.1), 



L(x,y) = 



Q h I (x)hu(y) 
x-y 

h I (y)h II (x) q 
x-y 



This matrix defines a bounded operator in £ 2 (X) if, for example, the operator T 
defined in (4.3) is £ 2 -bounded. 

Let us assume that the operator (1 + L) is invertible. This is automatically true 
if hi and hu are real valued, then L* = —L, and —1 cannot belong to the spectrum 
of L. 

As in §5, the operator L is integrable with N = 2. Indeed, 

L( X ,y)= hix)9 ' {y) + hix)9M , X ,yeX, 
x-y 



where 



h r (x), x G X/, f 0, x G Xi, 



0, x G X/7, t hu(x), x G Xi/, 

0, x G X/, f /i/(x), x G X/, 



hu(x), x G X/7, { 0, x G Xj-j-. 

The jump matrix w(x) = —f(x)g(x) t of the corresponding DRHP has the form 



w(x) = 




x G X/, 
x G Xu. 

It is readily seen that the DRHP is equivalent to the following conditions: 

• matrix elements ran and ra2i are holomorphic in C \ Xu; 

• matrix elements rai2 and ra22 are holomorphic in C \ Xj; 

• ran and ra.21 have simple poles at the points of X//, and for x G Xu 

Res u=x ran(w) = -hj(x)m 12 (x), 
Res u=x ra 2 i(w) = -h 2 I (x)m 2 2(x); 

• TO21 and TO22 have simple poles at the points of X/, and for x G X/ 

Res u=x rai 2 (w) = -/ijj(x)mii(a;), 
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• m(u) ~ 1 as u — > oo. 

Let us indicate how this setting of the problem can be obtained from the con- 
tinuous case. A continuous Riemann-Hilbert problem is equivalent to a system of 
integral equations. For the special type of Riemann-Hilbert problems described in 
§5 the system takes the form 

f /if(x)mn(x) f h 2 j(x)m 2 i(x) 

/ dx = m 12 {y), / dx = l + m 22 {y), 

Jsj x-y J^j x-y 

h n {x)m 12 {x) f hj I (x)m 22 (x) 



L 



En x-y 



dx = 1 + mn(y), / ^ii( x )m 22 (x) ^ _ m (\ 

Jy.u x-y 



For the first two equations y belongs to C \ Ej, and for the last two y belongs 
to C\E/j. 

A natural discrete analog of this system is the following one 

Ehj(x)mn(x) . . h 2 j(x)m 21 (x) 
= m 12 (y), > = 1 + 77122(2/), 
x — y ^-^ x — y 

xEXi xEXi 

Y ^ ( " )mi2(x) = 1 + mil (y), y ^ (x)m22(x) = m 21 (y), 
^— ' x — y ^-^ x — y 

xEXu xEXii 

which is equivalent to our DRHP. 

This analogy suggests that continuous RHPs can be obtained as limits of DRHPs 
when X "approximate" E, or, vice versa, DRHPs can be obtained as limits of RHPs 
when the contour E split into increasingly small parts passing through the points of 
X (the last observation is due to P. Deift). In §8 we will provide an explicit example 
of a continuous RHP and a DRHP such that the discrete problem converges to the 
continuous one in a certain limit. 

Proposition 4.3 in our special situation takes the following form, cf. §5. 

Proposition 6.1. Let m be a solution of the DRHP stated above. Then the matrix 
K = L(l + L)~ x has the form {with respect to the splitting X = Xi U Xu) 

K(x,y) = 

h I (x)h I (y)(-rn 11 (x)rn2i(y)+rn 2 i(x)rn 11 (y)) hj (x)h II (y) (mil (x)m 2 2 (y) — m 2 i (x)m 12 (y)) 

x—y x—y 

h II (x)h I (y)(m22(x)m 11 (y)-m 1 2(x)m 2 i(y)) hg {x)h I j (y) { — m 22 (x)m 1 2 (y)+m 12 (x)m 2 2 (y)) 

x-y x-y 

where the indeterminacies of type jj on the diagonal are removed by the L 'Hospital 
rule: 



K(x, x) 



hj(x)(—m' 11 (x)m 2 i(x) + m 21 (x)mn(x)), x G Xi, 
h 2 II {x){-m' 22 {x)mi 2 {x) + m' 12 (x)m 22 (x)) , x e Xu. 



Proof. A direct application of Proposition 4.3. Note that detm(C) = 1 by Lemma 
4.7, and 

( TO *)-i=f m22 ~ m2 A. □ 
\-m 12 run J 
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7. The discrete Bessel kernel 

Now we return to the problem stated in §2. As was explained in the previous 
section, the kernel K = L(l + L) _1 , where L is given by (2.1), can be expressed via 
the solution of DRHP of the special form discussed above with 



X = Z' = Z + - , Xi = Z' + = | - , - , . . . | , Xu = Z'_ 



1 3 

2' ~2 ! 



0% - § 
hi(x) = — — , h H (x) = 



r(*+£)' n-x+$ 

It will be more convenient for us to use the parameter r\ = \J~Q instead of 6>; then 

^ (x) = i>Ti)' hn{x) = n^Y)- (7 - 1} 

One way to extract the information about a solution of a RHP is to reduce the 
problem to a RHP with a constant (or not depending on a certain parameter) jump 
matrix. Then one can compare the solution with its derivatives with respect to the 
complex variable or with respect to a parameter, see, e.g., [KBI, Ch. XV]. We will 
employ this approach for our DRHP. 

The asymptotic condition of the DRHP requires that m(u) ~ 1 as u — > oo. 
Assume now that the next terms of the asymptotic expansion of the solution m{u) 
at infinity are given by the relation 

?HC 

Here a, (3, 7, 5 are some functions of r\. The asymptotics here and below is under- 
stood in the sense that \u\ — > 00 so that dist(w, X) is bounded from zero. 

Note that the DRHP has an obvious symmetry: since Xi = —Xu and hi(x) = 
hu(—x), the changes 

u ^_ u ( m n( u ) rn 12 (u)\ < f m 22 (-w) -m 2 i(-u)\ ^ ^ 



777,21 (u) m 2 2(u) J \-mi2(-u) mii(-u) 

do not affect the problem (the minus signs in the off-diagonal blocks appeared 
because of the relation Res M=x f(u) = — Res u= - X f(—u)). This immediately implies 
that a = —5, (3 = 7, so that 

™<«)=fi ?Vif2 L)+0(>). (7.3) 



1 J u \ P —ol 



V, 



Consider a new matrix 



n(u) — m{u) 



rf 
77"' 



Observe that n(u) is invertible for « G C\I, because m{u) is invertible by Lemma 

A T 
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If m(u) is the solution of the DRHP with /i's given by (7.1) and asymptotic 
behavior given by (7.3) then n(u) is the solution of the DRHP with 



and asymptotics 



"(«) = I I n °,)+l( a a l)+0(±))(«: »], ™ (7.5) 



1 y -u \/3 —a y \w 2 y y \ 77 

Note now that the jump matrix data of the new problem given by (7.4) do not 
depend on 77, so that the derivative ^ = d v n satisfies the same DRHP with possibly 
different asymptotics. By Lemma 4.5, the matrix (9 T? n)n _1 has no singularities in 
C. 

The asymptotics of (9 r? n)n~ 1 can be easily computed from (7.5): 



?i 2 



Note that ct has disappeared from the asymptotics. By the Liouville theorem we 
conclude that 

This relation already provides certain information about n. Indeed, it easily implies 
that the matrix elements of n satisfy second order linear differential equations in 77. 
However, the coefficients of this equation are still unknown — they are expressed 
in terms of the function f3{rf). 

The second trick which is commonly used in conventional RHPs in such a sit- 
uation is differentiation with respect to the variable, if the jump matrix does not 
depend on the variable. Unfortunately, we cannot reduce our DRHP to one with 
jump matrix not depending on u. However, this obstacle can be overcome in the 
following way. 

Introduce a new matrix 



p{u) = n(u) ^ r(M +^ 





1 

r(-u+£) 



If m(u) satisfies the original DRHP, it can be easily seen that p(u) is holomorphic 
in C. The residue condition of the DRHP and the asymptotics (7.5) in terms of 
p(u) take the form 



/;(>•) = ( i • )( ( ! I). .reX. (7.7) 
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rtOHll ')+<>(*)) *™ " I (7.8) 



1) u\P -a) \u 2 J J \ v ~ r- 

as w — > oo. 

Note that the relation (7.7) implies that the matrix p(x) is degenerate at the 
points x G X. 

Consider the matrix 



p(u) 



Pn(w) ^12^)^ _ ( pu(u + l) -pi2(u+l) 

P2l(u) P22(u) J \-p2l(u~l) P22(u-1) 



Clearly, p satisfies the condition (7.7), and (7.8) implies that, as u — > oo, 
pniu) = (1 + an" 1 + 0(n" 2 )) ; = - ' O^" 1 ), 

r(u+f) r(«+2) 
p 12 («) = + 0(u- 2 )) r l 1 ~ U ~\, = - f v ~]_ n ( § + o^" 1 ) 

= -(/to" 1 + 0(«- 2 )) = — ^ (-£ + Oiu- 1 )) 



2' 
,-tt+l 



p 22 («) = (1 - an" 1 + 0(n" 2 )) J 7 = 77 O^" 1 ). 

r(-u+f) r(-u+2) 

This means that the matrix 

satisfies the same DRHP as n(-u), and asymptotically 

n(it)n _1 (u) <■ 
By the Liouville theorem we get 



o £ 

_£ o 



u — > oo. 



^ , 

In the language of p and p this implies that 

P P 

Pll(w) =pn(w + 1) = -P2l(«), P2l(«) = -P2l(« - 1) = Pll(w), 

rj n 

P P 

Pl2(u) = -P\2{U + 1) = ~P22{u), P22(u) = P22(u - 1) = --pi 2 (lt), 
/ \ 2 

for all « 6 C. Hence, f^j =1, and /3 = ±rj. Then the equation (7.6) takes the 
form 

dr,n(v)=(r« T l)n(u). (7.10) 



(7.9) 
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Note that the change of signs of the off-diagonal elements of n swaps the solutions 
of (7.10) corresponding to (3 = r\ and (3 = —rj, so it is enough to consider the case 
(3 = —7]. Then (7.10) leads to the following relations 

u u 
d v nn(u) = -nn(it) + 2n 2 i(u), d v n 2 i(u) = --n 2 i(u) - 2nn(u). (7.11) 

This gives second order differential equations on nn('u), rai 2 (u): 

+ 4) n u («) = 0, (\ 2 - ^-±12+4) n 21 (u) = 0. 



u(u — 1) 



V v 2 



n 2 \\ 



General solutions of these equations are expressed through the Bessel functions 

Hi{u) = Vv (const 1 J u _i(2i]) + const 2 J_ u+ i(2i])^ , 
i( u ) = VV (const 3 J u+ i (277) + const 4 J_ u _i(2i])^ . 

The first line of (7.9) implies that (recall that (3 = —rj) 

const 1 = — const 3 , const 2 = — const 4 . 
The well-known differentiation formulas for the Bessel functions read 

d v (y/v J v (2rj)) =^-A j v (2r)) - 2^fq J„ +1 {2rj) 

- ' ' 2 J„{2i 1 ) + 2^J u - 1 {2i 1 ). 



Vv 

Using these formulas in the first equation of (7.11) we see that it is equivalent to 

4 const 2 yjrj J_ u _ 1 (2//) = 0. 

Hence, const 2 = and 

ran (it) = const ^frj J u _i(2ra), n 21 (u) = - const 0? J u+ i (2ra) 

with the same constant (depending on it). Since J u {2r]) ~ rf jY{y + 1) as 1/ — > 00, 
the asymptotic relation (7.5) implies that as u — > 00, 

nn(u) ~ T ^ v^J u _i(2r/), ra 2 i(u) ~ -r ( w + ^ 0) J «+i( 2 l)' 

Moreover, ran (it) and ra 2 i(w) have the same singularities as the corresponding ex- 
pressions in the right-hand sides above (simple poles at the points of £//). Hence, 

ran (it) = T (u + ^ JrjJ u _i(2-n), n 21 (u) = —T (u + ^ »/n J„ , 1 (2n). 
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Similar manipulations with 77-12 and ri2i show that 

ni2(u)=r(-u+^) VTjJ_ u+ i(2r,), n 22 (u) = T (-u + ^ ^J_ U _,(2 V ). 



Thus, 

It is easily verified that this matrix satisfies the original DRHP (the verification uses 
the symmetry relation J_ n = (— l) n J n , n£Z, for the Bessel functions). Now it is 
immediately seen that the kernel K(x,y) of Proposition 6.1 constructed from the 
matrix m above coincides with the discrete Bessel kernel (2.2), and the Theorem 
2.1 is proved. □ 

It is worth noting what happens if we choose f3 = r\ above. Then we get 

m(u)- /n( Ju -^ 2v) - J -»+hPv)\(v-"nu+l) \ 

m(U, -^U + i(2»7) J-u-^v) J { r,»T{-u+i))> 

and this matrix does not satisfy the DRHP. However, it satisfies the problem up 
to the change of sign in the residue condition. This change of sign is equivalent to 
the multiplication of hi and hu by \/— or to the change of sign of L. Thus, the 
kernel K constructed from the matrix m is connected with L by the relation 

^ T, 

K = 



L — 1 

8. Z-MEASURES 

It seems that the DRHP considered in §7 does not admit a scaling limit transition 
to a continuous RHP. The purpose of this section is to provide another DRHP for 
which such a limit exists. 

This DRHP and its scaling limit describe the so-called z-measures on partitions 
and the spectral decomposition of the generalized regular representations of the 
infinite symmetric group. We refer the reader to the papers [KOV], [BOll], [B012], 
[B013], where this material is thoroughly explained. 

We start with the description of the DRHP. As in §7, we take X = Z', X/ = Z^_, 
X11 = 7U_. The functions hi, hu are as follows (we add the superscript 'd' to the 
notation for the discrete problem and the superscript 'c' to the notation for the 
continuous problem) 



hj(x) = V 



hi Ax) = 



r(-x + 



Here z and z' are two complex parameters such that (z + k)(z' + k) > for all 
k G Z; £ G (0,1) is also a parameter; (a)k = T(a + k)/T(a) is the Pochhammer 
symbol. Note that the limit 
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brings us to the DRHP considered in §7. 

The continuous RHP in the notation of §5 has the form 

E = R\{0}, E/ = R + , E /7 = R_, 

UCI \ (zz')? £+£ _* 

hj(x) = — x 2 e 2 , 

y/T(z + l)T(z' + 1) 

h c n (x) = — (-X) 2 e2. 

y/T(-z + l)T(-z' + l) 

It is not difficult to see that 

lim hj([(l - Ox]) = hj(x), lim h d n ([(l - Ox]) = h c n (x). 

Here [y] denotes the integer which is closest to y. 

In such a situation it is natural to say that the DRHP approximates the contin- 
uous RHP as £ — > 1. 

The solution for the DRHP was obtained in [B012]: 

' /• ( •- « + |; ^ — : : 



m (u) 



F [i_ Z; i_ 2 ' ;u+ | ; ^ N ) 

_V 2 iZV p ( z z '- _ ?y + i. _i_ \ 



Here F(a,b;c;v) is the Gauss hypergeometric function. 

The solution for the continuous RHP was computed in different language in [Bl], 
[B2], see also [BOH]: 

m c (u) = 

u~ z+ 2 +1 et W z+Z ' +1 z - z > (u) \fzzJ (-u) z+Z 2 1 e~%W - z - z >-i z . z > (—u) 

^ 2 ' 2 i 2 ' 2 

-Vzz 7 e^W z+z >-i L s_(u) (-u)^2 e~%W - z - z >+i , _^{-u) 

2 ' 2 2 ' 2 

where W Ky/1 (v) is the Whittaker function. 

The correlation kernels K corresponding to m d (u) and m c (u) are called the 
hypergeometric kernel and the Whittaker kernel, respectively. 

The convergence of the DRHP to the continuous RHP is established immediately 
using the relation 

lim F ( a, b; x; 1 — — ) = y~ ^ e%W -a-b+i a -b (y). 

x^+oo \ y J 2 ' 2 



We have 



lim m d ((l - £)u) = m c (u), 



see [B012], [B013] for an explanation of the representation theoretic meaning of 
this limit transition. 

It is worth noting that for the discrete problem described in this section the rela- 
tion K = L(l + L) _1 holds for all values of parameters, and L and K are bounded 
operators in £ 2 (Z'). However, in the continuous case the operator L becomes un- 
bounded if \$l(z + z')\ > 1, and one should be careful to define (1 + The 
spectral analysis of the kernels L and K in the continuous case has been done in 

\rw~l „ — „i„„ fDAni 
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